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First-order field-induced spin reorientation transitions in multisublattice intermetallic compounds 
are considered within an anisotropic Heisenberg model. Unlike previous works, only leading-order 
£\j ■ anisotropy constants of the sublattices are included. The exchange interactions both between the 

rare-earth (RE) and transition metal (TM) sublattices and within the RE sublattice play an im- 
portant role in the model. The latter interaction results in the occurrence of metastable states. 
One of the minima corresponds to the non-collinear magnetic structure where the RE subsystem 
is divided into two sublattices, and their moments deviate in opposite directions from the TM mo- 
ment direction. The possible FOMP scenarios for R2Fei7, HoFenTi, and TbMneSng systems are 
discussed. 
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I. INTRODUCTION 



The phenomenon of FOMP (First Order Magnetization Process) was discovered during last years in many inter- 
metallic compounds. Here belong a number of practically important compounds of transition metals with rare-earths 
(RE) like R 2 Fei 4 B @,§, Nd 2 Coi 4 B @, R 2 Fei 7 @-|, RFenTi RMn 6 Sn 6 @, Sm x _ x Nd x Co 5 fHH, and also 



some ferrites |15| ]. In some directions, magnetization jumps with increasing external magnetic field are observed in 
these systems. Usually this effect is treated phenomenologically as a result of competition of magnetic anisotropy 
(MA) constants of different order @Q. Presence of large high-order constants enables one to obtain local min- 
ima of free energy and thereby first-order transitions. However, such an approach meets with a number of difficulties. 
i^i - Large values of if( n>1 ) required are not confirmed by crystal-field calculations (see Refs. |lP,|l8|). Besides that, high 
order K^ n ' are proportional to high powers of magnetization and should rapidly decrease with increasing temperature 
T. This effect should be stronger in the case where K^ n > come mainly from RE sublattice: the rare-earth moments 
can become rather small at T well below the Curie temperature Tq- The reason is that the exchange interaction Jrt, 
which acts on RE ions from the transition metal (TM) system, may be considerably smaller than the interaction Jtt 
which determines Tc JT^]. Thus, the local minima explaining FOMP should exist only at low temperatures. However, 
FOMP is observed as a rule in a wide temperature region [00. The dependences K^(T) fitted from magnetization 



curves turn out to be rather complicated and sometimes seem to be artificial. 

The data discussed lead to the conclusion that effects of many-sublattice structure are important in systems demon- 
strating FOMP |T^ ). It is known that in a number of intermetallic RE systems (Tm 2 Fei7, RFenTi, TbMn 6 Sn 6 ) 
temperature-induced orientational transitions are caused by competing anisotropy of RE and TM sublattices (in hex- 
aferrites BaZn 2 _ x Co x Fei6 2 7 fL5| |, competition of Fe and Co sublattices plays a similar role). It is natural to suppose 
that this factor is essential for the field-induced transitions too. Although the many-sublattice case was considered 
i | earlier p0[ |, these treatments were also based on the idea of competion of high-order constants. 

In the present work we propose a new mechanism of FOMP in a multisublattice magnet, which is connected with 
the formation of a non-collinear spin structure. The corresponding theoretical model is formulated in Sect. 2. In Sect. 3 
O ■ we consider the FOMP picture in various cases and discuss experimental situations. 



II. THE THEORETICAL MODEL 

We start from the Heisenberg model of uniaxial multisublattice crystal with inclusion of magnetic anisotropy. The 
MA constants in the problem are usually treated as phenomenological parameters which are obtained by fitting 
magnetization curves. When describing FOMP one introduces as a rule large high-order constants, ~ ~ 
To avoid the difficulties discussed above, we include only leading-order MA constants for different magnetic 
sublattices, Ki = Using the mean-field approximation we write down the expression for the total energy 

e = j2 J v cos (^ - 9 j o + Yj ( k * sin2 6i - HM *) ■ w 

i>j i 



1 



Here i,j are sublattice indices, Mj are the corresponding magnetic moments which make the angles 9i with the z 
axis, Jij are exchange parameters, H is the external magnetic field. For the sake of convenience, we have included 
the dimensionless factors of Mj = (Mj) into Jy, so that Jy oc Mj,Mj. The crystal and magnetic structures contain 
TM and RE sublattices which are coupled by the exchange interaction Jrt. The TM subsystem yields usually main 
contribution to magnetic moment, and the RE one to MA. The most strong interaction Jtt plays the role of a mean 
field in the TM subsystem and needs not to be included explicitly. 

As demonstrate crystal- field considerations, formation of non-collinear magnetic structures is possible in multisub- 
lattice structures, e.g., for the systems R2F914B The simplest case of a two-sublattice magnet was analyzed 

earlier in Ref. p9| . It was shown that, in the case where At and Ar have opposite signs (competing anisotropy), 
the sublattice moments become non-parallel provided that 

|Jrt(1/Ar + 1/A t )| <2 (2) 

If the absolute values of At and Ar are close, this condition can hold even for rather strong exchange Jrt. 

The two-sublattice problem with non-collinear magnetic structure enables one to obtain a local energy minimum 
with including magnetic field, but for one field direction only. To avoid this difficulty, we have to introduce, instead 
of one RE subsystem, two equivalent RE sublattices, RI and RII, with Ari = Aru = Ar, Jrit = Jriit = Jrt- 
This model preserves the correct axial symmetry and yields the energy minimum for both field directions. Of course, 
for real compounds with complicated crystal structures each 'sublattice' corresponds to several crystal positions. To 
stabilize the non-collinear structure with split RE sublattice we include the antiferromagnetic exchange interaction 
between the sublattices, Jrr > 0. The Jrt and Jrr values cannot be reliably obtained from existing experimental 
data, but for the long-range RKKY exchange interaction | Jrr| may be assumed not small. The essential property of 
such a model is the presence of metastable states already at H = 0. Owing to this fact, first-order phase transitions 
are possible which occur at crossing the energy levels in the external magnetic field. 



III. FOMP SCENARIOS FOR VARIOUS EXPERIMENTAL SITUATIONS 



The idea of the FOMP mechanism can be illustrated on a simple toy model where the TM moment is rigidly fixed by 
strong MA in the x direction (or in the xy plane for the three-dimensional problem) . The RE ions which possess more 
weak easy-axis MA tend to incline to the z direction, each ion having two local minima. The interaction Jrr results 
in that the degeneracy of the minima is lifted, so that the moments of RE sublattices deviate in opposite directions 
(Fig. la). With increasing magnetic field along the z axis, a first-order transition occurs to the configuration where 
RE moments are almost parallel (Fig. lb). A smooth transition does not occur because of the influence of MA of RE 
ions. 

Now we pass to a more realistic situation of finite (and even small) MA of TM subsystem. First we consider the 
situation where Ar > 0, A't < 0, and the easy-axis MA of RE sublattice predominates. This is the case for the 
systems Tni2Fei7 and TbMngSng, as follows from the existence of easy-axis - easy-plane spin-reorientation transitions 
with increasing temperature (|Ar| decreases more rapidly than |At|). Two possible metastable states are shown in 
Fig. 2, the exchange Jrt being supposed antiferromagnetic, which corresponds to heavy rare-earth ions, e.g., to Tb 
and Tm. In the case where Jrt < 0, which occurs for light rare-earths, e.g., for Nd and Pr, the picture differs by the 
opposite direction of RE moments. 

In the zero magnetic field the first state (Fig. 2a) is realized with the collinear spin configuration along the z axis 
(easy direction for the RE sublattice). At increasing the field H || x, the magnetization vector starts to rotate with 
retaining the nearly antiparallel orientation of TM and RE moments. For simplicity, we neglect in the analytical 
treatment the non-collinearity of TM and RE moments for this state, which is weak due to large value of | Jrt | ■ Then 
the energy in the state (a) is expressed in terms of the angle 9 beween the TM moment and z axis, 

A a = -2| Jr T | + Jrr - H{M T ± 2 Mr) sin0 
+ (A T + 2A R ) sin 2 (9 

= -2|J rt | + Jrr- (MT±2Mr) > (3) 
1 1 K 4(2A R + A T ) V ' 

the sign +(— ) corresponding to Jrt < (Jrt > 0). The second state (Fig. 2b) posseses a non-collinear spin con- 
figuration with split RE sublattice owing to the gain in the MA energy for the RE subsystem, and also due to the 
interaction Jrr. In this state the spin configuration depends weakly on magnetic field and the TM moment remains 
parallel to x axis, so that 
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£"b = — 2|Jrt| cosa + Jrr cos 2a 

-H(M T ± 2M R cos a) + K T + 2K R cos 2 a 

= K T -HM T -J RR - ^-Xy 2 (4) 

where a (0 < a < 7r/2) is the angle of the RE moment deviation from the x axis. The total energy in the state (b) 
has an almost linear field dependence and decreases more rapidly than that for the state (a) with rotating moment 
(Fig. 3). The angles 9 and a are obtained by minimization of Eqs.(||), (^), 

. H{M T ±2M R ) \J RT ~HM R \ 

Sm ° = 2(2K R + K T ) > C ° Sa = 2(J RR + K R ) > (5) 

Thus the local minimum (non-collinear state) occurs at 

J RR + K R >\J RT -HM R \/2. (6) 

Note that the minimum can occur with including magnetic field, even if it is a absent at H — 0. The transition field 
and the angle a (but not the magnetization jump) turn out to be small provided that 2\J RR + K R \ ~ | Jrt|- Since Jrr 
and K R are proportional to M 2 (T) (unlike the constants K^ n>1 ' in standard FOMP considerations, which decrease 
rapidly with increasing T) , the criterion (^|) depends rather weakly on T. Therefore the mechanism discussed enables 
one to obtain FOMP in a wide temperature region. 

The numerical calculations were performed by direct minimization of total energy (1) without approximations 
discussed above. The moment values were taken to be Mp = 4, M R \ = M R \\ — M R /2 — 1, which corresponds roughly 
to Pi2Fei7 compounds (e.g., for Tb2Fei7 the moment of Fe subsystem exceeds by two times that of RE subsystem, 
Mib = 18.6^b,Mfo = 36.6^b per formula unit Q). At some critical field E & becomes equal to Eb (Fig. 3), and 
the magnetization demonstrates a vertical jump (Fig. 4). It should be noted that the transition field is much smaller 
than the exchange interaction Jrr, unlike the situation of Fig. 1. After the jump, the magnetization increases rather 
slowly because of large value of Jrt- On the other hand, the magnetization depends very weakly on the field H || z 
because of the strong interaction Jrt- Owing to the deviation from the antiparallel alignment, the magnetization in 
the hard direction x after the transition turns out to be larger than the zero-field magnetization in the easy direction 
z (M x > M z = 4 — 2 = 2). Such a behavior is demonstrated by experimental dependences M(H) in a number of 
systems, e.g., in Tm 2 Fei 7 @,§. 

The zero- field canting angle a of the RE sublattice for the parameters of Fig. 3 makes up about 20°. Such small 
values do not contradict likely experimental data. Indications of formation of the non-collinear structure after the 
transition are given by neutron scattering data on Tni2Fei7 M and TbMngSng pi] ]. Non-collinear structures with the 
split RE sublattice were also observed for R 2 Fei4B systems p2| . 

In the system Tb2Fei7 the Fe sublattice has the easy-plane magnetic anisotropy, as well as the Tb sublattice 
(although Tb ions in different positions are supposed to have MA of different signs). Thus, the competing anisotropy 
situation does not occur. Nevertheless, in this case our approach also enables one to obtain FOMP, but the transition 
field turns out to be somewhat larger. Fig. 5 shows the result of calculation of magnetization for the constants of the 
same sign, K R > 0, > (the easy-plane and easy-axis problems are equivalent provided that we restirict ourselves 
to the plane model). One can see that before the jump the field dependence of the magnetization is practically linear, 
unlike Fig. 4. The state with split RE sublattices occurs only with increasing magnetic field. 

The curves obtained demonstrate a similarity to experimental ones. Let us try to perform a more detailed compar- 
ison with experimental data. The critical FOMP field is determined by the overall scale of the anisotropy parameters. 
Unfortunately, MA constants are usually fitted from magnetization curves with including high-order constants which 
contradicts our model. Therefore we shall concentrate on the cases where the parameters of the model can be de- 
termined from independent experiments. Using the data for Tni2Fei7 from Ref. Mrm — 12//b,Mf — 30/is 
per formula unit, JrmFe = 60J/cm 3 , Kp c = —2.5 J/cm 3 , i^Tm = 7 J/cm 3 , we can fit the experimental FOMP 
field, [j,qH = 5T, by taking Jrr/ Jrt — 0.5. Unfortunately, the value of the magnetization jump cannot be reliably 
determined from the results of 

The values of MA constants for Tb 2 Fei 7 can be estimated from NMR data ||| (see Ref. @), K Th = -10 J/cm 3 , 
Krb = 3 J/cm 3 for the positions 2b and 2d, -KTfc = — 3 J/cm 3 . Of course, the accuracy of these values is rather 
poor. After averaging the values of Xxb over sublattices such a relation between MA constants corresponds roughly 
to the parameters of Fig. 5. The rescaled FOMP critical field from this Figure turns out to be by about two times 
larger than the experimental value which makes up fioH ~ 4T. The agreement can be made better provided that 
we take somewhat smaller values of \K-p c \ and increase the ratio Jrr/ Jrt- Besides that, the situation in Tb2Fei7 is 
complicated because of existence of two FOMP-type transitions M. 
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Now we treat the case where easy-axis anisotropy Kt dominates over MA of RE sublattice, which has an opposite 
sign (such a stituation takes place, e.g., for the compound HoFenTi). Then the qualitative picture of FOMP turns 
out to be different from Fig. 2. As shown in Fig. 6a, at H = RE moments deviate from the axis z (easy direction for 
TM sublattice) by the angle 



arccos 



\Jm 



2( Jrr - K R ) 



(7) 



This configuration begins to rotate when including the small magnetic field in the x direction. At further increasing 
the field, the nearly collinear spin configuration shown in Fig. 6b becomes favorable. As follows from (Q), FOMP is 
possible at Jrr — 1£r > | Jrt|/2. Note that the solution (local minimum), corresponding to the non-collinear structure, 
can vanish with increasing H somewhat above the critical field of FOMP. A structure of the type Fig. 6a (but with 
ferromagnetic interaction Jrt < 0) was proposed for Yi-^Nd^Cos on the basis of magnetic measurements p5| . 

After the jump, the magnetization in the x direction continues to increase rather rapidly due to rotation of TM 
moments (Fig. 6b), and becomes saturated only when the moments lie along the x axis (M x = 2V. Absence of the 
saturation immediately after the jump corresponds to FOMP-II according to the classification [ p6p^ ] . The result of 
magnetization calculation is shown in Fig. 7. For simplicity, we restrict ourselves again to the plane problem (xz plane 
model). The parameter values are chosen to demonstrate existence of the saturation region, although FOMP can be 
obtained for smaller \Kij Jrt| and at smaller H too. The dependence obtained is in a qualitative agreement with that 
for HoFenTi |0]. This system is close to instability with respect to the formation of the 'easy-cone' structure |ll]] 
owing to competing MA of Ho and Fe sublattices with < 0, Kt > 0. 



IV. CONCLUSIONS 



To conclude, we have proposed a simple model which reproduces qualitatively the FOMP pictures in some RE 
systems with different sign combinations of sublattice MA constants. It should be noted that the field- induced 
transitions considered are similar in some respects to the spin- flop transitions in anisotropic antiferromagnets. In 
the latter situation, formation of metastable states and magnetization jumps are possible too (see, e.g., Ref. p6|). 
However, our picture is more complicated since this includes the TM sublattice. 

The consideration performed enables one to avoid some difficulties and reduce the number of fitting parameters in 
comparison with the treatments that use high-order MA constants. Further experimental investigations of magnetic 
structures, especially with the use of neutron scattering, would be of interest to verify the FOMP scenarios proposed. 

When considering the situation in other compounds, we meet with some problems. The description of FOMP for 
the systems Pr 2 Fei4B, Nd 2 Fei4B should also include non-collinear structures but this problem is more difficult. In 
the compound Nd2Fei4B a conical structure occurs in zero field, although the signs of MA for Nd and Fe sublattices 
coincide. Usually such a behavior is attributed to the influence of higher-order anisotropy constants; possibility of 
different role of various Nd positions is also discussed. It should be noted that the FOMP picture is similar for the 
compound Nd2Coi4B jj) where the Co ions have easy-plane MA. For the system Snii-^Nd^Cos Jl3| the RE sublattice 
already includes two subsystems with opposite MA signs. However, to obtain FOMP in spirit of our mechanism one 
still has to divide this sublattice into two symmetric parts. 

The author is grateful to Yu.P. Irkhin, N.V. Mushnikov, A.N. Pirogov, A.S. Ermolcnko and E.V. Rosenfeld for 
helpful discussions. The work is supported in part by the Grant of RFFI No. 00-15-96544 (Support of Scientific 
Schools). 



4 



Figure captions 

Fig.l. Schematic picture of FOMP in a simple toy model with ferromagnetic exchange Jrt, ^r > and large 
easy-plane Kt < (a) the spin configuration in the zero magnetic field (b) the spin configuration after the transition 
induced by the magnetic field in the x direction. 

Fig. 2. Schematic picture of FOMP in the case of antiferromagnetic exchange Jrt, -Kr > (a) the spin configuration 
in the zero field (b) the non-collinear spin configuration after the magnetization jump induced by the magnetic field 
H\\x. 

Fig. 3. The field dependence of the total energy in units of Jrt = 1. Other parameter values are = 0.1, Kt = 
—0.15, Jrr = 0.5. Long-dashed line corresponds to the state (a) in Fig. 2, and the short-dashed line to the state (b). 

Fig. 4. The field dependence of the magnetization in the x direction for the same parameters as in Fig. 3. The 
observable dependence demonstrating FOMP (vertical jump) is shown by the solid line. 

Fig. 5. The magnetization curve in the x direction for Jrt = 1, -Kr = 0.1, Kt = 0.1, Jrr = 0.6. 

Fig. 6. Schematic picture of FOMP in the case of antiferromagnetic exchange Jrt, i^R < 0. The spin configuration 
is non-collincar before the jump (a) and becomes collinear after the jump (b). 

Fig. 7. The magnetization curve in the x direction for Jrt = 1,Kr = — 0.2, Kt — 0.4, Jrr = 0.45. Long-dashed 
line corresponds to the state (a) in Fig. 6, and the short-dashed line to the state (b). 
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